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ABSTRACT 

A significant number of X-ray binaries are now known to exhibit long-term "su- 
perorbital" periodicities on timescales of ~ 10 - 100 days. Several physical mechanisms 
have been proposed that give rise to such periodicities, in particular warping and/or 
precession of the accretion disk. Recent theoretical work predicts the stability to disk 
warping of X-ray binaries as a function of the mass ratio, binary radius, viscosity and 
accretion efhciency, and here we examine the constraints that can be placed on such 
models by current observations. 

In paper I we used a dynamic power spectrum (DPS) analysis of long-term X- 
ray datasets (CGRO, RXTE), focusing on the remarkable, smooth variations in the 
superorbital period exhibited by SMC X-1. Here we use a similar DPS analysis to 
investigate the stability of the superorbital periodicities in the neutron star X-ray 
binaries Cyg X-2, LMC X-4 and Her X-1, and thereby confront stability predictions 
with observation. We find that the period and nature of superorbital variations in 
these sources is consistent with the predictions of warping theory. 

We also use a dynamic lightcurve analysis to examine the behaviour of Her X-1 
as it enters and leaves the 1999 Anomalous Low State (ALS). This reveals a signifi- 
cant phase shift some 15 cycles before the ALS, which indicates a change in the disk 
structure or profile leading into the ALS. 

Key words: stars - X-rays: binaries :pulsars raccretion disks rprecession periods 



1 INTRODUCTION 

In a small but growing number of bright X-ray binaries a 
third (or "superorbital") period is present in addition to 
the usual periodicities. Broadly speaking, superorbital vari- 
ations divide into two observational classes, with clear, sta- 
ble X-ray intensity variations in the ~30-day range forming 
one category (e.g. SMC X-1, Her X-1, LMC X-4), and more 
quasi-periodic, longer periodicities in the ~ 50-200 day range 
forming the other (see Clarkson et al, 2002, hereafter paper 
I, and references therein). Modulations in the longer-period 
class are distinctly "quasi-periodic" in that even long-term 
monitoring datasets do not yield precise periods in their 
power spectra, but a broad peak, often superposed on an 
apparent "red-noise" continuum. Examples of these include 
Cyg X-2 (P ~70-80d; Paul, Kitamoto & Makino 2000, here- 
after PKMOO), GX354-0 (P ~70d; Kong, Charles & Ku- 
ulkers 1998), X1916-053 (P ~83d; Homer et al. 2001) and 
X1820-30 (P ~171d; Chou & Grindlay 2001). 

With the exception of X1820-30, it has become custom- 
ary to interpret superorbital periodicities from both classes 



as due to a precessing, possibly warped accretion disk, with 
diflterences in amplitude due to diflterences in inclination. 
This scenario dates from the early history of X-ray astron- 
omy, in which a warped disk was proposed to explain the 
behaviour of Her X-1 (Gerend & Boynton 1976, Petterson 
1977). The intense radiation pressure from the central X-ray 
source drives the disk warp, and tidal forces lead to preces- 
sion on the observed 35-day cycle. The precession causes 
periodic obscuration of the central source, giving rise to a 
modulation of the X-ray flux at the precession period. This 
model has been the subject of intense theoretical examina- 
tion in recent years (see e.g. Wijers and Pringle 1999, or for 
the current state-of-the-art, Ogilvie and Dubus 2001, here- 
after ODOl, and references therein). 

In Paper I it was argued that varying absorption could 
not account for the superorbital periodicity in SMC X-1, as 
variation in the 1.3 - 12.1 keV band was mirrored at higher 
energies. Instead, the warp could manifest itself in two ways: 
as a variation in uncovered emitting area (as in Gerend and 
Boynton 1976), or as a variation in accretion rate onto the 
neutron star. 



2 Clarkson, Charles Coe, Laycock 



A further manifestation of disk warping not mentioned 
in Paper I is modification of disk throughput by X-ray irrar 
diation of the disk. X-ray irradiation from the central source 
influences the disk temperature and viscosity (van Paradijs 
1996, King & Ritter 1998). In a flat disk, self-screening 
strongly diminishes this effect (Dubus et al 1999), but if the 
disk is at all warped, irradiation reaches beyond the inner 
disk in an aaimuthally asymmetric way. In this case irradi- 
ation can significantly heat the disk (Hynes et al 2002) and 
is thought to affect the rate of mass flow through the disk 
(Dubus et al 2003), with measurable time- varying effects on 
the disk temperature profile and central source intensity. 

Scenarios other than disk warping arc possible, how- 
ever. We merely introduce these mechanisms here and refer 
the reader to section 6.2 of Paper I for more information:- 
(i) Also postulated to explain the superorbital periodicity 
in Her X-1 was the possibility of precession of the magnetic 
axis of the neutron star (Triimper et al 1986; Lamb et al 
1975). Here, the variation of accretion from equator to pole 
causes a torque that leads to the precession of the mag- 
netic axis, (ii) A third body can induce an eccentricity in 
the XRB orbit, causing orbital precession and decay. The 
precession manifests itself as a change in accretion rate on 
a period determined by the orbital period of the XRB and 
the third body (Chou & Grindlay 2001). This process has 
been postulated to explain the ~ 171 day period in X1820- 
303 but for SMC X-1 the low amplitude of modulation of 
the pulse arrival times (Wojdowski et al 1998) requires an 
extremely finely tuned triple system, so this mechanism was 
discarded for that system (paper I). Mechanisms for which 
no prototypical system have yet been identified involve (iii) 
variations in the extent to which the donor overfiows its 
Roche Lobe as a result of stellar pulsation (see, e.g: Welilau 
et al 1992), or (iv) varying location of an accretion X-ray 
bright spot (as happens in the UV in CV's, e.g: Rutten, van 
Paradijs & Tinbergen 1992). 



X-ray binaries, the majority of the torque exerted by the 
mass input is close to rc, as even in LMXB with large accre- 
tion disks, the matter stream is believed to continue to the 
vicinity of the circularisation radius. Accordingly, we follow 
ODOl in taking the location of mass addition to be ~ rc- 
Figure 7 shows the division of XRB into regions of stabil- 
ity to warping from this analysis. The location of SMC X-1 
at the boundary between regions stable to simple mode 
warping and combinations of modes supports our assertion 
that the long-term X-ray behaviour of SMC X-1 is due to a 
competition of warping modes. 

The location of the boundaries between stability regions 
in rt - q space depends strongly on the global disk viscosity 
parameter, a and the neutron star accretion efficiency rj. 
Uncertainty in these values translates into uncertainty as to 
the location of the boundary between regions and thus the 
stability classification of a system. However, the adopted 
values of a = 0.3 and = 0.1 are usually assumed to be 
appropriate for neutron star XRB (Frank, King & Raine 
1995, ODOl and references therein). 

In this paper we examine the long-term behaviour of 
a selection of X-ray persistent neutron star XRB in which 
superorbital periodicities are well-established. We apply the 
techniques from paper I, in particular the DPS, to explore 
the phenomenology of these variations. We focus on neutron 
star XRB for which persistent superorbital periods have pre- 
viously been reported. Our goal is to relate the long-term 
behaviour of the sources to their position in the - q dia- 
gram, so as to directly confront the stability predictions of 
the radiation-driven warping framework with observation. 
If succesful, it will establish the time evolution of periodic- 
ities as a powerful indicator of XRB properties in systems 
whose properties are less well-understood, for example the 
Ultracompact X-ray Binaries (Clarkson et al 2003, in prep) 



1.1 Accretion Disk Warping 

In Paper I we presented the techniques of examining super- 
orbital variations in XRB, focusing on SMC X-1, the source 
which shows the clearest variation in its superorbital pe- 
riod. The dynamic power spectrum (or DPS) was used to 
chart the behaviour of the superorbital period. An appar- 
ently cyclic variation (on an even longer timescale of ^7 
years) in the superorbital period was evident in the DPS. It 
was suggested that the behaviour of SMC X-1 is consistent 
only with a precessing, warped accretion disk, where a com- 
petition of warping modes leads to variation in the warp, 
giving rise to the shifting superorbital period. 

This interpretation gains support from the theoretical 
stability analysis of ODOl, which extends a-disk theory to 
radiation-driven warping, with correct analytic expressions 
for the torque on a disk element. Solutions are examined nu- 
merically to determine stability to warping for initially flat 
disks, steadily precessing disks and nonlinear warps. They 
find that for reasonable values of the global disk viscosity 
a and accretion efficiency q, the stability of an accretion 
disk to warping is determined by the binary separation rt 
and mass ratio q = A further variable is the location of 
the mass input to the disk, whether close to the Lindblad 
resonance point or the circularisation radius Vc- For most 



2 RXTE OBSERVATIONS 

Launched in December 1995, the Rossi X-Ray Timing Ex- 
plorer (RXTE) carries an All Sky Monitor (ASM), which 
gives regular coverage of the entire sky. Typically 5-10 read- 
ings - called "dwells" - are taken of each of a list of sources 
per day, lasting about 90 seconds per dwell. Timing informa- 
tion is provided to within a thousandth of a day, as well as 
crude spectral information. The ASM is sensitive to photon 
energies between 1.3 and 12.1 keV, broken into three energy 
channels (1.3-3.0 keV, 3.0-5.0 keV and 5.0-12.1 keV). Nearly 
7 years of data (January 1996 - August 2002) from the ASM 
were used in our analysis. For all analyses reported in this 
paper, the ASM data were selected by background and qual- 
ity of coverage; all points with background level above 10 
cs~^ were rejected as part of a filtering technnique modeled 
after Levine et al (2000). Crude X-ray spectral information 
can be obtained from the ASM data by constructing hard- 
ness ratios of the three channels. In this work we define the 
ratio to be the sum of the count rate in the two high energy 
channels (3.0 - 12.1 keV) to that at low energy (1.3 - 3.0 
keV). 
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3 SOURCE SELECTION 

Each source used in this analysis was selected on the basis 
of three criteria. Firstly, independent estimates of r/, and 
q must be available, which restricts our search to optically 
identified systems with full orbital solutions. Secondly, the 
sources must be persistent X-ray emitters at a high enough 
level to provide long-term datasets of sufficient quality to 
detect long-term modulations, which favours neutron star 
XRB. Third, each system used must show independent evi- 
dence for a persistent accretion disk. 

Table 1 shows the sources selected for this study, to- 
gether with the warping behaviour predicted by ODOl. The 
strong dependence of the stability boundaries on a and r) 
leads us to choose systems with similar compact objects in 
the hope that for such systems these parameters will take 
similar values. Poor coverage of the - q diagram forces us 
to select systems with both high and low - mass companions. 

Figure 1 shows sections of the RXTE/ASM Ughtcurves 
from the sources used here, binned appropriately to im- 
prove signal-to-noise. The liglitcurve of Cyg X-2 shows sim- 
ilar morphology to that of SMC X-1, though at count rates 
about a factor of 20 higher. Its morphology is also more 
complex than a single periodic cycle, which leads to the de- 
tection of more than one periodicity (PKMOO). The ASM 
lightcurve for Her X-1 shows a clear variation on the well- 
known 35 day cycle, with a secondary peak half a cycle after 
the main period of extended activity. LMC X-4 also shows 
evidence of a cycle that is superficially similar to that of 
SMC X-1, although the lightcurve shows more scatter. 



3.1 Discarded Sources 

Several well-constrained, persistent neutron star XRB are 
within the regions of interest in the Vh — q diagram, but we 
choose not to analyse them here as they arc unlikely to show 
significant long-term cyclic phenomena in the RXTE/ASM. 
Circinus X-1 is well within the region of the stability diagram 
corresponding to mode 1 and higher warping, and is a lumi- 
nous neutron star XRB (Tennant 1986). However, its orbit 
is strongly suspected to be highly eccentric, among many 
other properties the system shares with A0538-66 (Charles 
et al. 1983). As pointed out in ODOl, the stability analysis 
as it stands is not valid for non-circular orbits. 

4U 1907-1-09 is a pulsar XRB with massive compan- 
ion (e.g. van Kerkwijk et al 1989) in a wide orbit, that has 
shown a monotonic increase in pulse period since its dis- 
covery (Mukerjee et al 2001). Transient QPO's have twice 
been discovered at frequencies (~10 - 400 mHz), the inter- 
pretation of which requires an accretion disk, however small. 
This phenomenon suggests that temporary formation of an 
accretion disk is the mechanism responsible for the spin- 
down, and that this must occur during the majority of or- 
bital cycles (Mukerjee et al 2001). The interval over which 
this disk persists is small compared to the 8.4-day orbital pe- 
riod (in 't Zand et al 1998), suggesting that radiation-driven 
warps would not persist long enough to be measurable by 
RXTE/ASM. 

X2127+119 has recently been shown to consist of two 
sources which cannot be resolved by RXTE (White & An- 
gelini 2001; Charles, Clarkson & van Zyl 2002). While re- 
solving the mystery of how a source can show evidence for 



both an ADC and Type I bursts, it is not possible to relate 
components in the observed RXTE/ASM X-ray fiux to a 
single system. 

OA01657-415 shows spin-up and spin-down trends 
without any correlation with X-ray luminosity > 20keV 
(Baykal 2000), suggesting no persistent accretion disk is 
present. In the case of OA01657-415, these episodes can- 
not be explained by wind accretion (Baykal 1997), but are 
thought to arise through temporary formation and destruc- 
tion of a disk, in the manner implied by Bildstcn et al (1997). 
A small correlation between accretion rate and pulse fre- 
quency change in the RXTE waveband is taken as evidence 
that this is indeed the case. 



4 DYNAMIC POWER SPECTRUM 

Data were analysed using our DPS approach, designed to be 
sensitive to long-term quasi-periodic variations. The Lomb- 
Scarglc (LS) periodogram code (Scargle 1982, 1989) was 
used in conjunction with a sliding 'data window' to pro- 
duce power density spectra (PDS) for a series of overlap- 
ping stretches of the time series. Adjustable parameters in 
the analysis were the length of the data window and the 
amount of time by which the window was shifted each time 
to obtain the overlapping stretches of data. The code ac- 
counted for variations in the number of datapoints per in- 
terval such that the power spectrum resolution was identical 
for each interval. 

Choice of data window length was infiuenced by the 
limits of reliable period search with the PDS; we choose 
two complete cycles as the minimum needed for a reliable 
detection. The range of claimed supcrorbital periodicities 
runs from ~ 30 - 200 days, so we set our window length 
at 400 days. We choose the maximum length that meets 
our criteria because this brings about the sharpest period 
detections. 

4.1 Statistical Significance of Detected Periods 

Once a period is detected, its statistical significance must be 
addressed. For most of the sources in this paper, long-term 
variability is superposed on frequency-independent white- 
noise. The statistical significance of a detected periodicity 
over this noise is determined by computing the LS power 
spectra of 10^ simulated 400-day datasets (assuming gaus- 
sian white-noise) with the sampling of the true lightcurve, 
then measuring the fraction of datasets that produce the LS 
power of interest or greater (figure 2). For the 400-d data 
windows used in this analysis, an LS power of 9 corresponds 
to a 99.9% confidence level that the detected periodicity is 
real, and not a collusion of random noise and data sam- 
pling. These sources show LS powers in excess of 30, sub- 
stantially in excess of even the 99.9% significance level. The 
highly coherent nature of the variation profiles (figure 1) ar- 
gues against frequency-dependent red-noise (e.g. Timmcr & 
Konig 1995, Homer et al 2001), affording confidence in our 
significance levels. 

The error in the detected period is determined empiri- 
cally. A sinusoidal variation at the detected periodicity, with 
amplitude corresponding to the detected peak LS power is 
superimposed on white noise, with random sampling on the 
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Figure 1. Sections of the RXTE/ASM lightcurves of the sources examined in this work, binned appropriately to emphasise their 
long-term behaviour. 
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Figure 2. Significance and period tests for 400-day RXTE/ASM datascts. Left: Probability that a maximum power LSmax can be 
generated by white- noise alone as a function of LSmax- Right: Number of trials producing detected period P. The period error is the 
3a width of the fit to the distribution combined with the magnitude of any shift between the most likely period detection and the input 
periodicity. 
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Table 1. Neutron star XRB exhibiting superorbital modulations 



Name 




1 


P 1 

sup 


Mx 


^orh 


-M2 


Prediction ^ 




(xlO«GMi/c2) 




(d) 


{A^0) 


(d) 


(xlO-*MQ2/r-i) 




Her X-1 


3.1 


1.56 


35 


1.36 


1.700 


1.4 


Persistent mode 


LMC X-4 


4.5 


10.6 


30 


1.35 


1.408 


3.5 


Persistent mode 


SMC X-1 


8.9 


11.0 


50 


1.6 


3.890 


2.3 


Mode dominant but other modes possible 


Cyg X-2 


8.0 


0.34 


78 


1.78 


9.844 


2 


Complex combination of modes 0,1 and higher 



Claimed superorbital period ^ Based on position in figure 7 of ODOl ^ See text for references 



RXTE/ASM count rate and over the same length of time 
as a data window from the DPS. The LS periodogram of 
this dataset is performed and the peak frequency measured. 
This process is repeated for 10® simulated datasets, and the 
resulting spread of detected periodicities about the input 
periodicity measured. The detection accuracy is quoted in 
units of 3cr of the fit to the resulting distribution in period 
detections. In one case (Cyg X-2) the peak period detection 
was shifted from the input period by an amount correspond- 
ing to the beat between the detected period and ~ 8 times 
the time interval of a data window, which has been incorpo- 
rated into the quoted error on the period detection (table 2 
and figure 2). 



5 RESULTS 
5.1 Her X-1 

Her X-1 is the prototypical disk warping system, with a 
'^35-day periodicity detected at optical and X-ray wave- 
lengths. This periodicity is almost certainly due to a pro- 
cessing, warped accretion disk, as first suggested by Gerend 

6 Boynton (1976) to explain the optical lightcurve, and cor- 
roborated by the evolution of the pulse profile (Scott et al 
2000), and variations in X-ray illumination of the compan- 
ion HZ Her (Leahy & Marshall 1999, Leahy, Marshall & 
Scott 2000). The accretion disk is thought to be tilted and 
twisted (Gerend & Boynton 1976, Schandl & Meyer 1994, 
Schandl 1996). With X-ray pulsations at 1.24s, eclipses at 
1.7 days and an A7 Roche lobe-filling companion (Middled- 
itch & Nelson 1976) of ~2 A/©, the mass ratio and binary 
radius are constrained (but see Leahy & Scott 1988 for un- 
certainties due to possible subsynchronous donor rotation). 
On three occasions, the source has failed to reach a high 
state at the time expected from the 35-day cycle (Parmar et 
al 1985; Vrtilek et al 1994; Oosterbroek et al 2000). These 
Anomalous Low States (ALS) last for several 35-day cycles 
each, and are thought to be the result of a state change in 
the accretion disk (Still et al 2001). 

The DPS of Her X-1 (figure 3) shows an extremely sig- 
nificant (>>99.99 %) periodicity at 34.98 days, along with a 
less significant secondary peak at 17.5 days. This secondary 
peak is not an artefact of the search technique; as figure 3 
shows, the DPS has picked up the secondary peak at phase 
0.3 in the 35-day cycle. Her X-1 is the only system suspected 
of disk precession to show this secondary feature so signif- 
icantly. Her X-1 also shows the clearest relation between 
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Figure 3. DPS of Her X-1. Contours spaced at intervals of LS 
power 100. 
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Figure 4. Superorbital lightcurve (left) and count rate-hardness 
relation (right) for Her X-1. The spectral hardness varies with 
the ASM count rate and there is a clear relation between the two 
quantities. 



ASM count rate and spectral hardness of all the sources 
studied here (figure 4) . 



5.2 LMC X-4 

LMC X-4 is a ~13.5 s X-ray pulsar with a 15 M© 0-type 
companion that eclipses the X-ray source every 1.4 days. 
QPO's have been found at 2-20 mHz during large X-ray 
flares (Moon & Eikenberry 2001) confirming the presence 



6 Clarksoriy Charles Coe, Laycock 



Period (Days) 
70 60 50 40 



Period (Days) Period (Days) 

Q 100 80 70 60 50 40 30 25 UK) 8(1 70 60 50 41) 30 25 




E 

Q 



0.01 0.015 0.02 0.025 0.03 0.035 

Frequency ( 1 / day) 

Figure 5. DPS of LMC X-4. Contours spaced in intervals of 
LS power 10 (by comparison, LS power 9 corresponds to 99.9% 
significance over noise - section 4.1) 
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Figure 6. Superorbital lightcurvc (left) and count rate - hardness 
relation (right) for LMC X-4. Though a more noisy dataset, there 
is a correlation between the variations of the two quantities. 



of orbiting matter at leeist ~ 0.25 rdrc from the neutron 
star. Disk modeling of the orbital and processional optical 
lightcurvc (Hccrnskcrk & van Paradijs 1989) suggest the ac- 
cretion disk may extend out to about 4rc. 

The DPS of LMC X-4 (figure 5) shows the simplest be- 
haviour of the sources analysed here: a periodicity at 30.28 
days that is both stable and persistent throughout the en- 
tire dataset. The DPS is not sensitive to the ~ 0.02% yr~^ 
systematic variation of this periodicity reported recently by 
Paul & Kitamoto (2002). At LS power ~ 38, this 30-d cy- 
cle is somewhat less significant than for Her X-1 (LS Power 
~ 300), but the power still corresponds to > 99.9% sig- 
nificance. There appears to be slight jitter in the peak pe- 
riod, but at ± ~ 0.05 days this is not significant when com- 
pared with the 3a uncertainty for this source of ±0.46 days. 
LMC X-4 shows a relationship between RXTE/ASM spec- 
tral hardness and count rate (figure 6) that is somewhat less 
steep than for Her X-1. 

5.3 Cygnus X-2 

Although not a pulsar, Cyg X-2 nevertheless has a well- 
constrained binary orbit, with a 0.5 - 0.7 Mq companion in 
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Figure 7. DPS of Cyg X-2. Contours spaced at LS powers of 200. 



a 9.84-day period (Casaxes, Charles & Kuulkers 1998), un- 
usually long for a LMXB. Type I X-ray bursts (e.g: Kahn 
& Grindlay 1984) established Cyg X-2 as a neutron star 
XRB. Modehng of lUE spectra (Vrtilek et al 1990) yields a 
best-fit outer accretion disk radius of ~ 0.6 Rl (the Roche 
lobe radius) , while modeling of optical lightcurves (Orosz & 
Kuulkers 1999) suggests the disk radius may be closer to 
0.9 Rl- a pattern of high- and low- intensity X-ray states 
has been seen in early RXTE/ASM, Vela 5B and Ariel-5 
observations (Wijnands et al 1996), on timescales of ~77.7, 
~77.4 and ~ 77.7 days respectively, while the Ginga/ASM 
lightcurvc shows the strongest periodicities on periods of 
53.7 and 61.3 days (PKMOO). This variation has been at- 
tributed to radiation-driven accretion disk precession (Vr- 
tilek et al. 1997). 

Cyg X-2 shows complex behaviour in the DPS, with 
an extremely significant (LS > 900) periodicity at period 
~ 59 days, seen over the period from day number 600-900 
(figure 7). This feature seems to form as two separate pe- 
riodicities converge, then diverge as the feature vanishes. 
There are also significant recurring features at ~ 40 days. 
The ~ 77-day periodicity noted previously may be present 
during the beginning of the RXTE/ASM dataset, though is 
of equal or lower significance than its harmonic. In any case 
the longterm variation of Cyg X-2 is not the reliable clock 
suggested by Wijnands et al (1996). In contrast to Her X-1 
and LMC X-4, the binned superorbital lightcurvc of Cyg X-2 
shows a negative relation between the RXTE/ASM spectral 
hardness and count rate (figure 8), as expected if the su- 
perorbital lightcurve is due to an obscuration effect (section 
6.3). 



6 DISCUSSION 

6.1 Her X-1: Inclination Change 

The clear hardness-intensity relation in the RXTE/ASM 
lightcurve (figure 4) is not by itself an indication of the 
physical mechanism (Tanaka 1997). In fact we find no fea^ 
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Her X-1: Dynamic Lightcurv 




Figure 8. Supcrorbital lightcurvc (left) and count rate - hard- 
ness relation {right) for Cyg X-2, over day numbers 600-900. In 
contrast to Her X-1 and LMC X-4, the spectral hardness appears 
to have a negative correlation with the count rate. 
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Table 2. Detected stable periodicities 



Name 


Psup (d) 


AP,„p (d)i 


Her X-1 


35.0 


0.16 


LMC X-4 


30.3 


0.46 


Cyg X-2 


59.0 


0.45 



3(T 



ture of the lightcurve of this source that contradicts the 
precessing, warped accretion disk scenario of Gerend and 
Boynton (1976). As the warp precesses, it could uncover 
regions deeper towards the central accretor, predicting the 
correspondence between spectral hardness and count rate 
observed. The asymmetry in the supcrorbital lightcurvc can 
be brought about by a combination of the system inclination 
and vortical extent of the warp (Gcrcnd & Boynton 1976, 
Schandl & Meyer 1994, Schandl 1996). 

The superorbital periodicity is completely absent from 
a significant fraction of the data, marking the third ALS de- 
tected in this source to date (Oosterbroek et al. 2001). This 
is not an effect of data sampling, as might be expected for 
a source with such a low count rate; on the contrary, the 
persistent emission at ~ 1 cs~^ is still well sampled. Mea- 
surements by BcppoSAX (Oosterbroek ot al 2001) show that 
the most recent ALS is accompanied by a rapid spin-down of 
the neutron star, in which the spin-down rate is some nine 
times faster than the spin-up rate during "normal" X-ray 
output. At the same time, optical observations during the 
ALS show (Margon et al 1999) that the companion is still 
strongly irradiated. Spectral fits and lightcurve modeling of 
RXTE/PCA pointings during the ALS suggest the X-ray 
output during the ALS consists of X-ray reflection from the 
companion (Still et al 2001). The current interpretation for 
this behaviour is that the line of sight to the central source 
has been obscured during the ALS, causing the turn-off, by 
an increase in line-of-sight column density due to the disk. 

This would impact the supcrorbital lightcurve in two 
ways. Firstly, the asymmetry between the peaks in the fold 
would be altered, but whether as an increase or decrease 
would depend on the current inclination and the precise disk 
profile. Secondly, a phase shift in the superorbital fold would 
be seen during the change of inclination. The twist of the 
disk brought about by the warp means that any given feature 
in the warp would be seen at different disk azimuth angles 
for different heights above the nominal flat disk equator. 



Figure 9. The 35-day dynamic lightcurve of Her X-1. The phas- 
ing of the peaks undergoes a shift before the onset of the Anoma- 
lous Low State (ALS), and the relative significance of the peaks 
alters. When the ALS ends, the phasing is unaltered from the 
early years of the dataset. 



In figure 9 we plot the dynamic lightcurve (DL) of Her 
X-1. In a manner analogous to the DPS, this breaks the 
dataset into windows of sufficient length to give good cov- 
erage, then folds these windows on the ephemeris specified. 
Because fewer cycles per window are needed for detection in 
a fold than a power spectrum, the windows can be shorter 
for the DL than for the DPS; it was therefore decided to use 
independent windows. 

From the DL we see that there is indeed a phase shift 
taking place over some fifteen cycles before the ALS. The 
secondary peak in the 35-day cycle appears to decrease in 
relative amplitude just before onset of the ALS. The cut-off 
in the superorbital period detection, however, is abrupt, tak- 
ing place over three cycles at most. By the end of the ASM 
dataset, the system shows behaviour similar to the begin- 
ning of the dataset, and the secondary feature is recovered. 

This phase shift and profile change is most naturally in- 
terpreted as a change in the line-of-sight disk profile. There 
are two good candidates for the nature of this change. It has 
been suggested (Wijers & Pringle 1999, van Kerkwijk et al 
1998) that if a disk is sufficiently warped, the inner disk may 
actually flip over and start to precess in the opposite direc- 
tion to the outer disk. This would produce the extremely 
rapid spin-down observed during the ALS as the direction 
of flow of accreting matter would then oppose the spin of the 
neutron star. Following this interpretation, wo would expect 
to be observing the neutron star through the accretion disk 
for at least some of the ALS, which would produce a harder 
spectrum for the ALS than for the normal state. However, 
the continuum X-ray emission is not spectrally harder for 
the ALS than for the emission prior to and after the ALS 
(flgure 10). A further difficulty is that one would expect to 
see some manifestation of the change of disk orientation in 
the overall ASM lightcurve. 

A second interpretation involves the disk changing its 
inclination with respect to the observer. We suppose here 
that the X-ray intensity from a twisted, warped disk vaxies 
as the amount of emitting area uncovered in the inner disk 
region varies. If such a disk were to change its inclination 



8 Clarkson, Charles Coe, Laycock 



ON: Count Rote vs Hardness 




LfvIC X-4: Dvnamic Lightcurv 



ASM c/s 

ALS: Count Rate vs Horaness 




Figure 10. Her X-1 Hardness-Intensity relation for the normal 
state (top) and the ALS {bottom), using 1-day data bins to im- 
prove S/N. There is no significant difference between the low- level 
emission in the two states. 



with respect to the observer, the uncovered area at each 
point in the cycle would also change with the inclination. 
A twisted disk would naturally produce a variation in the 
phase of features in the profile. That the ALS is accompanied 
by rapid spin-down of the neutron stax presents a difficulty 
for this scenario, unless either (i) the accretion flow at the 
magnotosphoro/inner disk boundary changes direction to act 
against the neutron star rotation, or (ii) the spin of the 
neutron star causes ejection of matter in the manner of e.g: 
4U1907+09 (section 3.1), producing a rapid spin-down as 
the accretion flow lessens due to a reduction in accretion 
rate through the disk. In scenario (ii), the spin-down during 
the ALS would suggest that the Her X-1 disk is reducing 
its inclination to the binary plane. The extent to which the 
uncovered area is sensitive to the inclination change, and 
thus the level of other observable effects of this change, is not 
at present clear. A further effect might bo the variation with 
the warp of the area open to X-ray heating, which might 
bring about a reduction in mass throughput depending on 
the timescale for such heating (c.f. Dubus et al 1999). 



6.2 LMC X-4: Stable variation 

At some five times more distant than most of the sources 
in this study, LMC X-4 exhibits low count rate and hence 
lower S/N, as can be seen in its 31-day dynamic lightcurve 
(figure 11), which furthermore undergoes no significant shift 
in phase or overall intensity over the ASM dataset. The su- 
perorbital lightcurve itself is somewhat different in shape to 
that of Her X-1 (figures 4 and 6), with the primary peak 
in the profile less symmetric, and no significant secondary 
peak is detected. A secondary peak might be present, how- 
ever: assuming a similar relative amplitude as that in Her 
X-1, such a peak would occur at only ^0.2 c/s. 

This suggests that if the accretion disk in this source 
does shift its inclination, as suggested for Her X-1, the 
LMC X-4 accretion disk either does not do so at all on the 
timescale detectable with the DL, or the disk profile and 
inclination are such that any existing wobble does not im- 
pact the superorbital fold to the same extent. We can place a 
rough upper limit on superorbital phase variation of features 
in the dynamic lightcurve of ~ ±0.1 orbit. 





Figure 11. Dynamic lightcurve of LMC X-4. Data window length 
chosen to include 5 cycles per window - each bin in this plot is 
then the average of typically 30 datapoints. 



6.3 Cyg X-2: Multimode Variation 

For Cyg X-2, in which the DPS behaviour does not show 
a single stable periodicity, such analysis is not possible. 
Indeed, its behaviour is highly complex, with significant 
periodicities forming and shifting over the 7 years of the 
ASM dataset. This explains the range of modulations re- 
ported for this source over the years (e.g: Smale & Lochner 
1992, PKMOO). Of particular interest is the apparent con- 
vergence, reinforcement and subsequent divergence of peri- 
odicities about April 1996 - January 1999 (figure 7). Because 
the accretion disk is a fluid entity and the DPS overlap be- 
tween data windows of the same order as the warp preces- 
sion periods, we do not expect too discontinuous a change 
in the variation shape: in other words we expect the pattern 
(primary peak-|-single harmonic) to persist over a timescale 
of several data windows. At the point at which the strong 
DPS peak appears to split into two modes, we interpret the 
mode evolving towards lower frequencies as the dominant 
frequency, as it is more significant. It is thus easy to in- 
terpret this feature as a single periodicity with ill-defined 
phasing (hence the broad peaks in the DPS), with a period 
that changes in a manner superficially similar to SMC X-1, 
only on a much shorter period: ~2 years compared to the 
~7 years of SMC X-1. 

However, the behaviour may be more complex. Closer 
examination of the secondary peak at ~ 40 days suggests it 
is not a simple harmonic of the primary period, as its period 
is some 20% longer than the expected location of the har- 
monic (see also PKMOO) . Furthermore, this structure is still 
present over the last 1500 days of the dataset, even though 
the fundamental it would be a harmonic of has all but dis- 
appeared. This feature may represent either (i) a slightly 
shifted harmonic that represents a deviation from a purely 
sinusoidal variation profile, or (ii) a mode of variation inde- 
pendent of the ~ 59-day variation. 

We can test which scenario describes Cyg X-2 better 
through simulations of varying periocicities. We take the 
superorbital lightcurve over day number 600-900 (from Jan 
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Figure 12. DPS of simulated Cyg X-2 lightcurves. Top: input 
data consists of the superorbital lightcurve from day number 600- 
900, which undergoes period change in a manner approximating 
the fundamental from figure 8. Middle: As for top, but with 40- 
day periodicity added that is slaved to the fundamental, to ap- 
proximate a non-sinusoidal shape of the variation. Bottom: 40-day 
periodicity kept steady as the fundamental varies. 
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Figure 13. Cyg X-2: DPS of the Ariel-5 dataset, binned to one- 
day averages, with contours at LS powers of 9 (corresponding to 
99.9% significance over white-noise). 



6 1996), then scale the period over time in a manner ap- 
proximating our identified dominant periodicity. The DPS 
of this simulation is quite diflForont in character to the real 
DPS (figure 12); the most obvious differences arc its far lower 
significance and frequency than the ~ 40-day periodicity ac- 
tually seen. A 40-day periodicity was then artificially added 
but slaved to the fundamental to simulate a non-sinusoidal 
variation. The DPS structure that results varies more closely 
with the dominant periodicity than that actually seen. If the 
~ 40 day periodicity is made completely independent of the 
fundamental, however, by being held steady as the funda- 
mental varies), the DPS that results is similar in morphology 
to that actually seen (figure 7). 

We thus suggest here that the apparently variable pe- 
riodicity at ~ 59 days and the more stable periodicity at ~ 
40 days represent independent variations, rather than a sin- 
gle variation of non-sinusoidal shape. This suggestion is in 
agreement with the conclusion of an earlier analysis based on 
the static periodogram of the early years of the RXTE/ASM 
(PKMOO). This suggestion gains some support when we no- 
tice that the 40-day periodicity is close to the expected har- 
monic of the ~ 77-day periodicity reported in the archival 
Vela 5B and Ariel-5 datasets. To attempt to constrain the 
properties of this behaviour over a longer baseline, the Ariel- 
5 and Vela-5B datasets were subjected to DPS analysis. The 
Vela-5B dataset was used in the filtered, binned form pro- 
vided by NASA HEASARC. Enough points were present to 
use the static periodogram to confirm the 77.4-day period 
already noted elsewhere (e.g: Smale & Lochner 1992), but 
not enough to follow its time evolution in a statistically sig- 
nificant way. However the Ariol-5 dataset has higher S/N, 
allowing less coarse binning and thus providing enough dat- 
apoints to permit DPS analysis. As figure 13 shows, the peri- 
odicities are highly variable over time, producing the highly 
complex structures seen in static periodograms (PKMOO). 
The ^ 77-day periodicity is strongly detected over the mid- 
dle 400 days of the dataset, while its harmonic at ~ 40 
days is detected throughout most of the ASM dataset. How- 
ever, as with the RXTE/ASM dataset, there is a further 
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periodicity which appears variable and transient. Thus the 
apparent combination of two or more separate variations 
does not appear to be a new behaviour. Furthermore, this 
complex behaviour has now been confirmed using two inde- 
pendent datascts separated by 20 years. 

The X-ray manifestation of the precessing warp may be 
different from that evidenced from the canonical disk pre- 
cessing systems SMC X-1, Her X-1 and LMC X-4. Cyg X-2 
shows a decrease in spectral hardness with count rate (figure 
8), which at first sight would be inconsistent with a process- 
ing warp. However, the variation of uncovered emitting area 
invoked to explain the superorbital periods in Her X-1 and 
LMC X-4 requires the column density through which the 
bulk of the X-rays roach RXTE to remain the same as the 
warp processes. If, however, the emitting area is permanently 
obstructed by the outer accretion disk, assumed not to be a 
significant emitter in the ASM bandpass, the variation can 
be brought about by changing the column density through 
which the majority of the X-rays travel. In this case, the 
warped disk acts to subtract flux from the central source, 
rather than contributing significantly to any emission. This 
removal of flux will be seen as X-ray obscuration, which 
thereby causes the observed spectrum to harden. 

6.4 Stability of Accretion Disks to Warping 

We now use our results for the other systems showing su- 
perorbital modulations to test the stability framework of 
ODOl. We first attempt to identify the warping modes cor- 
responding to the superorbital variations identified in this 
work. With their stable behaviour in the DPS and simple 
superorbital lightcurvcs. Her X-1 and LMC X-4 show clear 
evidence for mode warping. SMC X-1 shows a simple su- 
perorbital fold that varies over time, probably due to disk 
instabilities brought about as other modes begin to form 
(paper I). Cyg X-2 shows at least two independent modes of 
warping, thought to represent mode 0,1 and possibly higher 
warp modes. 

The sources thus behave in ways qualitatively predicted 
by the stability analysis of ODOl (table 1). Her X-1 and 
LMC X-4 are both predicted to show stable mode warping, 
which is indeed observed. SMC X-1 is thought to show warp- 
ing commensurate with its position near the mode and 
mode 1 regions. Cyg X-2 is deep within the region predict- 
ing mode 1 and higher modes, which is consistent with the 
appearance and interaction of at least two warping modes. 

We illustrate the evolution of disk warping with binary 
radius in a schematic bifurcation diagram, shown in figure 
14. This charts the evolution of procession solutions as the 
control parameter is increased. According to the ODOl 
predictions, for low values of the control parameter rb, no 
radiation-driven disk warping is predicted at present. In the 
region corresponding to stable mode predicted warps, such 
warps arc indeed found, with the precession frequency in- 
creasing with rb (Her X-1 and LMC X-4 type variations). 
Near the boundary between regions, modes begin to com- 
pete, producing an instability in the superorbital period 
(SMC X-1 type variations). Finally, once the boundary into 
full multi-mode variations is crossed, strong periodicities 
form and interact (Cyg X-2 type variations). When we over- 
plot the location of these sources on the schematic, along 
with their DPS results, we find that the long-term behaviour 



of the sources agrees qualitatively with the predictions of ta- 
ble 1, as can be determined from the DPS of the sources. 



7 CONCLUSION 

We have analysed the time variation of superorbital pe- 
riodicities present in three persistent neutron star XRB: 
Her X-1, LMC X-4 and Cyg X-2. Her X-1 shows a shift 
in phasing of the superorbital period shortly before enter- 
ing its third Anomalous Low State, which we interpret as 
evidence for a change in disk inclination. LMC X-4 shows 
stable mode warping, with little change over the period 
of the RXTE/ASM dataset. Cyg X-2 shows a combination 
of at least two separate modes of variation, suggesting more 
complex warping shape than for the other sources. Further- 
more this behaviour has persisted for at least ~ 1000 binary 
orbits, as shown by new analysis of the Ariel-5 dataset. 

We have identified the most likely warping modes to 
produce the observed variations in Her X-1, LMC X-4 and 
Cyg X-2, and together with SMC X-1 have constructed an 
observational stability sequence based on their behaviour in 
the DPS. We find that this qualitatively agrees with the 
predictions of the stability analysis of ODOl and establishes 
the DPS as a powerful tool to probe XRB properties. 
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